Amplification of Local Instabilities in a Bose-Einstein Condensate with Attractive 

Interactions 
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We study the collapse of large homogeneous Bose-Einstein condensates due to intrinsic attractive 
interactions. We observe the amplification of a local instability by seeding a momentum state p 
and suddenly switching the scattering length negative via a Feshbach resonance. As required by 
momentum conservation, we also observe the appearance of atoms in the conjugate momentum 
state. The time scale for this depletion process is found to be comparable to that for global collapse, 
implying that this process will be the primary decay channel for large homogeneous condensates. 
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Our current understanding of the collapse of Bose- 
Einstein condensates (EEC's) with attractive interac- 
tions is incomplete. While the experiments in ^Li pro- 
vided many insights into the formation kinetics and sta- 
bility , it was not until the discovery of externally in- 
duced Feshbach resonances that it became possible 
to tune the value of the scattering length and study in 
detail the effects of an attractive mean field potential. 
This technique was used by Donley et al to study *^Rb 
with a negative scattering length, and they have observed 
rich dynamics inherent in the collapse of the conden- 
sate. Among the most intriguing observations was the 
formation of low energy 'bursts' and 'jets' which were 
ejected out of the condensate While theoretically the 
enhancement of quantum fluctuations could give rise to 
such phenomena |Q, ^, ^ , there is currently no consensus 
on the exact mechanism by which it occurs. 

So far, theoretical developments have been limited by 
having only one experimental testing ground. The ex- 
periments using ®^Rb were done using small condensates 
15000 atoms), where the attractive mean field energy 
H is comparable to or less than the fuv level spacing of the 
harmonic trapping potential. In this paper, we study the 
collapse of large sodium condensates far in the Thomas 
Fermi regime >> fiw), where the spatial profile of 
the condensate is relatively homogeneous. Much of the 
dynamics of such a system is then described by local 
phenomena. When the interactions become attractive, 
Yurovsky ||^ predicts that local instabilities with mo- 
mentum on the order of the (imaginary) speed of sound 
will undergo exponential growth. Simultaneously, mo- 
mentum conservation requires atoms to be generated in 
conjugate momentum states. Since amplification hap- 
pens on the time scale of the chemical potential h/ ^, the 
resulting quantum evaporation of the zero momentum 
condensate atoms can happen faster than the global col- 
lapse for large condensates, where the whole condensate 
'implodes'. 



We probe this decay channel by seeding a particular 
momentum state with an initial population, then sud- 
denly switching the scattering length negative via a Fes- 
hbach resonance. At the same time, the trapping poten- 
tial is turned off so all subsequent dynamics are due only 
to the intrinsic attractive interactions. The resulting am- 
plification and the associated generation of atoms in the 
conjugate momentum state verifies the theory. We end 
with a discussion on the different collapse time scales of 
competing processes and show that for large condensates, 
this decay channel becomes dominant. 

The theoretical basis for the amplification of local in- 
stabilities is the dispersion relation for the elementary 
excitations in a Bose-Einstein condensate: 
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where hq is the density, U — Airh^a/m = /i/no is the 
contact potential, a is the scattering length and m is the 
mass. For an elementary excitation whose momentum 
satisfies |pp/2r7i < 2|/i|, an instability forms when /i < 
(i.e. a < 0) and oscillatory behavior gives way to expo- 
nential growth or decay. A formal derivation gives the 
evolution of these low momentum modes as 
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where A(p) is given by |e(p)|/;i (Eq0) when a < § 
and is the destruction operator for mode p. The in- 
stability of the mode pair (p, — p) results in correlated 
growth, where the creation of an atom in the -l-p state 
is accompanied by the creation of an atom in the — p 
state. A similar phenomenon is also responsible for the 
four wave mixing process observed in 0, Hl- At higher 
momentum, the energy becomes real again as the excita- 
tions now have enough kinetic energy to stabilize them 
against the attractive interactions. 

In our experiments, we created large cigar-shaped 
sodium condensates in the F ~ 1, mp — —1 spin state 
with typical atom numbers of ^^30 million and peak den- 
sities of 3 X 10^^ cm"-^ in a loffe-Pritchard (IP) magnetic 
trap. Following this, they were adiabatically loaded into 
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a mode-matched 1064 nm cylindrical optical dipole trap 
with trapping frequencies of 250 Hz by 2 Hz and held 
for 1 s to allow transient excitations to damp out. In 
the next 117 ms, the axial magnetic field was ramped up 
to an intermediate value just short of the 1195 G Fes- 
hbach resonance 1^. In order to obtain a homogeneous 
field across the condensates, we used two pairs of coils 
coaxial with the condensates. The first pair (bias coils) 
provided a large bias field ('^1200 G) with some gradient 
(~6 G/cm). The second pair (pinch coils) produced a 
cancelling gradient together with a small bias field (~8 
G), resulting in a nominally uniform bias field. The scal- 
ing of the field with current was characterized at low 
fields, allowing us to adjust our fields to within 0.4 G. At 
this point, the condensates were slightly compressed by a 
small axial magnetic field curvature, giving them radial 
and axial dimensions of 15 /im and 1.5 mm respectively 
and a speed of sound c = ^JUn^lm of 8 mm/s. 

We imprinted a low momentum excitation onto a sta- 
ble condensate (a > 0) using a two photon optical Bragg 
transition to couple atoms from the zero momentum state 
to a low momentum state p . The Bragg beams were 
pulsed on for 400 /is and were directed at an angle of 15° 
and 20° respectively from the long axis of the conden- 
sate, creating phonons propagating with a momentum of 
mx2.5 mm/s in a predominantly radial direction. With 
this momentum p, the minimum value of |a| required for 
amplification to happen is estimated to be —0.06 nm. 
The beams were red-detuned from the sodium D2 line 
by 3 nm to minimize Rayleigh scattering, and had a fre- 
quency difference of Aw — 2-k'k 700 Hz. 

At this point, we changed the scattering length a al- 
most instantaneously by switching off the small pinch 
coils, thereby stepping up the magnetic field by 8 G and 
entering the Feshbach resonance. This gave us a well- 
defined initial condition similar to that in We chose 
to use the 1195 G Feshbach resonance as it is the broadest 
resonance in sodium with a width of ~ 4 G and the scat- 
tering length crosses zero at a magnetic field below the 
singularity. This allowed us to tune the value of a contin- 
uously from the background scattering length oq to zero 
to arbitrary negative values by varying the intermediate 
field value at the end of the slow initial ramp. Since the 
pinch coils are much smaller than the bias coils, their mu- 
tual inductance is small and the sudden switch off only 
affected the bias field adversely by < 1 G. Moreover, the 
sudden jump reduced the time available for three-body 
decay, thereby overcoming the problematic high loss rate 
first observed by Inouye et al ||] . We also turned off the 
optical trapping potential in the same instant, allowing 
the system to evolve under its intrinsic attractive inter- 
actions. After a variable time of to 1 ms, the magnetic 
fields were turned off, causing the scattering length to 
become positive again. A pump pulse from F = 1 to 
F' = 2 was then applied and radial absorption images on 
the F = 2 to f = 3 cycling transition were taken after 
6 ms of ballistic expansion. These images provided the 
momentum distribution of the condensate. 
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FIG. 1: Amplification of unstable excitations in a BEC with 
a < 0. In (a), the unamplified +p excitations are shown mov- 
ing out of the condensate after 6 ms ballistic expansion. In 
(b), the -l-p excitations have been amplified and are accompa- 
nied by the formation of — p excitations after 600 of hold 
time at a = —0.82 nm. The field of view is 960 yLxn by 100 



The images in Fig. ^ probe the radial dynamics of the 
condensates and provide dramatic visual verification of 
quantum evaporation. Fig. |^(a) shows the -fp excitations 
moving out of the condensate without any amplification. 
In contrast. Fig. ^b) was taken after the condensate had 
been held at a = —0.82 nm for 600 ^s. Not only was the 
number of atoms in the -|-p momentum state significantly 
amplified, it was accompanied by the formation of excita- 
tions in the — p momentum state, here seen moving out of 
the condensate in the opposite direction. These observa- 
tions clearly demonstrate the instability of a condensate 
with negative scattering length. 

Due to the large aspect ratio of our condensates, only 
part of the condensate could be imaged at high mag- 
nification. However, since the -fp and — p excitations 
were created predominantly in the radial direction, this 
was not a limitation. The small 'kinks' that were also 
apparent in our condensate are most likely a result of 
imperfections in our optical dipole trap. Yet rather than 
degrade the signal, they highlight the parallel contours 
between the condensate and the ridge of excitations as 
atoms move out with a definite momentum. 

In order to perform more quantitative tests of this phe- 
nomenon, we first characterized the negative scattering 
length dependence on the field by directly probing the 
strength of the attractive interactions. To do this, we 
prepared optically trapped condensates close to a Fesh- 
bach resonance as above. The confining infrared laser 
beam was then replaced with a repulsive 3 nm blue- 
detuned "antitrap" beam as we simultaneously jumped 
to negative scattering lengths. At the right intensity, 
the antitrap beam provided the correct amount of repul- 
sive dipole force needed to compensate for the attractive 
interactions within the condensate and suppressed any 
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FIG. 2: Characterization of the Feshbach resonance. The 
light intensity needed to balance out the interatomic interac- 
tions gives the dependence of the scattering length on mag- 
netic field. The width of the resonance was found to be 
2.4±0.4 G. 

global contraction of size. However, this is an unstable 
equilibrium and any sloshing of the condensate or mis- 
alignment of the laser beam caused the condensate to 
be repelled. Therefore the antitrap was fine-aligned to 
milliradian accuracy such that a condensate with a > 
was ripped apart radially into a hollow cylinder. Using 
this method, we were able to stabilize an attractive con- 
densate significantly for 0.2 to 2 ms, depending on a{B), 
before unavoidable losses became significant. The radial 
dimension was used to monitor the mechanical dynamics 
of the condensate occurring on the 250 Hz time scale of 
the trap frequency. For minimal distortion of the spatial 
image, absorption images were taken after only 2 ms of 
ballistic expansion necessary for the high magnetic fields 
to die out. 

In equilibrium, a cc Fattractivc = ^Icpuisive cc I. There- 
fore we plot the light intensity / needed for stabiliza- 
tion vs. magnetic field B in Fig. || and obtain the scat- 
tering length dependence on the field. A red-detuned 
laser beam with a similar detuning of 3 nm was em- 
ployed to obtain the points in the positive scattering 
length regime. By fitting the expected Feshbach curve 
a{B) = ao(l -I- ) to our data, we find the width 

AB of the 1195 G Feshbach resonance to be (2.4±0.4) 
G. Here, uq = 3.3 nm is the triplet scattering length at 
high fields |l^ . While the exact position Bq of the res- 
onance is unimportant for our experiment, the width of 
the resonance is crucial as it determines the range of the 
magnetic field we have to work within ||ll|| . 

A quantitative analysis of the growth in the -|-p and 
— p modes was performed by monitoring their occupa- 
tion number as a function of hold time in the attractive 
regime (Fig. |^). Taking into account the high loss of 
atoms during this process due to three-body decay or 
inelastic two-body collisions, we normalize the number 
count for each mode to the number of atoms in the con- 
densate at the end of the hold time. The maximum du- 
ration of amplification was limited by the lifetime of the 
condensate, which was about 600 fis for a = —1.35 nm. 
Following Eq. 0, an exponential dependence was fitted to 



FIG. 3: Amplification of excitations in a BEC with a < 0. 
The growth of the normalized mode population in both +p 
and — p modes as a function of hold time at a = —1.35 nm is 
fitted to an exponential dependence according to Eq. ^ which 
gave a common growth rate of 5.89±0.83 ms~^. A variable 
offset was allowed to account for the initial seeding of the +p 
mode. 



the data, yielding a common growth rate of 5.89±0.83 
ms~^ for both modes. This agrees well with the theoret- 
ical growth rate of 5.57 ms~^, estimated using our initial 
mean field oi h x 5 kHz. 

We also investigated the dependence of the growth 
rate on the magnitude of the negative scattering length 
\a\. By varying the scattering length and extracting the 
growth rate of the excitations as above, we observed a 
strong increase of the growth rate with \a\ (Fig. ^ ). A 
fit to the theoretical prediction A(p) = &i\/|ap-62 (Eq. 
|l|) yielded the fit parameters h = (1.11 ± 0.14) x 10® 
nm^/^ms^^ and 62 = 0.078 ±0.04 nm, in agreement with 
theoretical estimates of 1.55 xlO® nm^/^ms^^ and 0.06 
nm respectively. For this sequence of measurements, the 
initial mean field was fj, = h x A kHz. The large error 
bars reflect the high sensitivity of the dynamics to the 
magnetic field. In particular, the three-body recombina- 
tion and inelastic loss rates as a function of a have not 
yet been well characterized, which limits the accuracy of 
our data. 

The results presented here prove conclusively that 
quantum evaporation is a part of the complex dynam- 
ics that take place during the collapse of an attractive 
condensate. While we select a particular mode for obser- 
vation, the effect is predicted to happen for all modes sat- 
isfying the condition |pp/2r7i < 2|/i|. At short times, the 
pairwise emission of atoms also implies that the number 
of atoms in the conjugate mode pairs will be exactly cor- 
related, although as the condensate becomes increasingly 
depleted, higher order effects will degrade the correlation 

i- 

Since quantum evaporation is intrinsic to the conden- 
sate, a natural question to ask would be how large a part 
it plays in the unperturbed collapse of the condensate. 
We study this by comparing the observed time taken for 
the condensate to decay completely without initial seed- 
ing (Fig. ^), with the theoretical prediction for the 
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FIG. 4: Growth rate of excitations in a BEC with a < 0. The 
rate increases with the magnitude of the negative scattering 
length \a\. The sohd hne shows a best fit to the predicted 
bi y^\a\ — 62 dependence. 



In addition, we compare the observed collapse time to 
the decay time predicted for global collapse. Even for 
condensates initially far in the Thomas-Fermi regime, 
their spatial profile is not completely uniform and the 
inherent pressure gradient will cause the condensate to 
collapse inwards. As the condensate compresses, the den- 
sity will increase and sharply enhance the three-body re- 
combination loss rate, which goes like Ug. We model 
the radial evolution using the root mean square radius 
i? = J r'^\ip{r)\'^rdr. For our cylindrical condensates 
with an aspect ratio of 100:1, an analytical solution 
for the resulting 2D dynamics exists ||l^, |l^, given by 
R = 4(-^ — uj^R) where E is the total energy of the sys- 
tem. Since ui — and a is constant throughout in our 
experiments, E is conserved and obtained from initial 
conditions. The time taken to reach i? = is 



o 
o 




^dccay — 



(4) 



12 3 4 

Scattering Length (nm) 



FIG. 5: Time scales for BEC decay. The observed collapse 
time is compared to the time taken for the condensate density 
to fall to 1/e of its original value (ti/e) due to instability 
amplification (solid curve, Eq. ^ and the expected idecay for 
global collapse (dashed curve, Eq. n). 
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From Eq. ^, we extract the time ti/^ taken for the con- 
densate density hq to decay to 1/e of its original value 
with our parameters (Fig. ||, solid line). As both the 
number of unstable modes and the amplification rate in- 
crease with a, there is a significant decrease of ii/g with 
\a\. We also considered the effect of the initial quantum 



depletion (~ 1%, 
impact is small. 
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Eqs. H and |4| are plotted in Fig. || and their intersec- 
tion separates the graph into two domains. In the first, 
a is small and global collapse is predicted to dominate 
over quantum evaporation. For higher a's, the converse 
is true as amplification rate and the number of unstable 
modes increase. The comparison with the data suggests 
that the decay of the condensate is not dominated by 
three-body decay, but is caused by amplification of un- 
stable modes. We were unable to study the condensate 
lifetime at even more negative scattering lengths since 
the condensate then decayed almost instantaneously. In 
order to further separate the two time scales, the global 
collapse time needs to be much slower compared to ti/e. 
This would require a combination of a weak trap with a 
high number of atoms, which is currently out of reach. 

In conclusion, we have shown that large condensates 
far in the Thomas Fermi regime undergo amplification of 
local instabilities when their scattering length becomes 
negative. We have studied the dependence of amplifi- 
cation rate on the magnitude of the negative scattering 
length and found reasonable agreement with the theory. 
For our parameters, this quantum evaporation process 
becomes comparable or faster than the global collapse. 
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